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ABSTRACT 

The spectral signatures of asymmetry in Type la Supernova (SN la) explosions are investigated, 
using a sample of late-time nebular spectra. First, a kinematical model is constructed for SN la 
2003hv, which can account for the main features in its optical, Near-Infrared (NIR), and Mid-Infrared 
(Mid-IR) late-time spectra. It is found that an asymmetric off-center model can explain the observed 
characteristics of SN 2003hv. This model includes a relatively high density, Fe-rich region which 
displays a large velocity off-set, and a relatively low density, extended 56 Ni-rich region which is more 
spherically distributed. The high density region consists of the inner stable Fe-Ni region and outer 
56 Ni-rich region. Such a distribution may be the result of a delayed-detonation explosion, in which 
the first deflagration produces the global asymmetry in the innermost ejecta, while the subsequent 
detonation can lead to the bulk spherical symmetry. This configuration, if viewed from the direction 
of the off-set, can consistently explain the blueshift in some of the emission lines and virtually no 
observed shift in other lines in SN 2003hv. For this model, we then explore the effects of different 
viewing angles and the implications for SNc la in general. The model predicts that a variation of the 
central wavelength, depending on the viewing angle, should be seen in some lines (e.g., [Ni II] A7378), 
while the strongest lines (e.g., [Fe III] blend at ~ 4700A) will not show this effect. By examining optical 
nebular spectra of 12 SNe la, we have found that such a variation indeed exists. We suggest that 
the global asymmetry in the innermost ejecta, as likely imprint of the deflagration flame propagation, 
is a generic feature of SNe la. It is also shown that various forbidden lines in the NIR and Mid-IR 
regimes provide strong diagnostics to further constrain the explosion geometry and thus the explosion 
mechanism. 

Subject headings: radiative transfer - nuclear reactions, nucleosynthesis, abundances - supernovae: 
individual (SN 2003hv) - supernovae: general 



1. INTRODUCTION 

Type la supernovae (SNe la) provide a powerful tool 
to investigate the cosmological parameters and the na- 
ture of the dark energy (Riess et al. 1998; Perlmutter 
et al. 1999). Thanks to the uniformity of their optical 
peak luminosity, once a phenomenological relation be- 
tween the light curve shape and the peak luminosity is 
applied, SNe la can be used as reliable cosmological dis- 
tance indicators (Phillips 1993; Phillips et al. 1999). 

There is a general consensus that SNe la are ther- 
monuclear explosions of a carbon-oxygen white dwarf 
(WD) (e.g., Nomoto et al. 1994, 1997; Wheeler et al. 
1995; Branch 1998; Hillebrandt & Niemeyer 2000). A 
Chandrasekhar-mass WD has been favored as a progen- 
itor for the majority of SNe la (e.g., Hoflich & Khokhlov 
1996; Nugent et al. 1997; Mazzali et al. 2007). 
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The thermonuclear runaway starts with the ignition of 
deflagration bubbles (e.g., Nomoto et al. 1976, 1984). 
It has been suggested that the deflagration flame may 
turn into a detonation wave (delayed detonation model, 
or deflagration-detonation-transition model; Khokhlov 
1991; Yamaoka et al. 1992; Woosley & Weaver 1994; 
Iwamoto et al. 1999; Ropke & Niemeyer 2007b), al- 
though the details of the transition have not been clari- 
fied yet. 

The first deflagration phase may well proceed in a very 
asymmetric way (Niemeyer et al. 1996; Garci'a-Senz & 
Bravo 2005; Livne et al. 2005; Jordan et al. 2008). The 
deflagration wave propagates under the work of the buoy- 
ancy force, and thus a small perturbation in the progen- 
itor structure could result in a global asymmetry. Rota- 
tion and convection in the progenitor WD could provide 
the seed for this asymmetric deflagration propagation. 
For example, there is a possibility that the convection 
in the progenitor WD is dominated by a dipole mode 
(Woosley et al. 2004), which likely results in highly off- 
axis ignition and propagation of the deflagration flame 
(Ropke et al. 2007c; Kasen et al. 2009). 

It is, however, observationally challenging to put con- 
straints on the geometry of the explosion. Measurements 
of polarization suggest that a large global asymmetry 
does not exist in SNe la (Wang et al. 1996), with only 
a few exceptions (Howell et al. 2001). However, the 
polarization probes mainly the outer regions of the ex- 
panding SN ejecta, at least with the existing telescopes 
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TABLE 1 

SN Ia Optical Nebular Spectra Sample 21 



SN 




Phase (Days) 


[Ni II] 


References 


1986G 


+ 1 03 

TWO, 


+257 


Yes 


Oistiani pt al | , 1992"l 


1990N 


+186, 


+227, +255, +280, +333 


Yes 


Gomez ct al. (1996) 


1991bg 


+1 1 7 


+199 


No 


Tiirattn pt al (IQQfi^ 

_1_ U.1 OjLi lr»_> <^ Lr Oil.. \±!y<j\Ji 


1991T 


+258, 


+284, +316 


No 


Gomez ct al. (1996) 


1992G 


+106, 


+128 


No 


Gomez et al. (1996) 


1994D 


+106 




Yes 


Gomez et al. (1996) 


1996X 


+246 




No 


Salvo et al. (2001) 


1 998an 


+211 


+231 +241 


Yes 


Branch ct al. (2003) 


1998bu 


+249^ 


+329' 


Yes 


Cappcllaro et al. (2001) 


2000cx 


+125, 


+147, +360 


Yes 


Li et al. (2001); Sollerman ct al. (2004) 


2001V 


+106 




Yes 


Matheson ct al. (2008) 


2001cl 


+398 




Yes 


Mattila et al. (2005) 


2002dj 


+275 




Yes 


Pignata et al. (2008) 


2002cr 


+216 




No 


Kotak ct al. (2005) 


2003cg 


+385 




No 


Elias-Rosa et al. (2006) 


2003du 


+138, 


+141, +209, +221, +272, +377 


Yes 


Anupama et al. (2005); Stanishev et al. (2007) 


2003hv 


+110, 


+143, +320, +358, +398 c 


Yes 


Leloudas et al. (2009); Gerardy et al. (2007); Motohara et al. (2006) 


2004dt 


+152 




Yes 


Altavilla ct al. (2007) 


2004eo 


+228 




No 


Pastorello ct al. (2007) 


2005cf 


+267 




No 


Leonard (2007) 



a The data are compiled from the SUSPECT data base, except for SNc 2000cx, 2001el and 2003hv. For SN 2003hv, NIR data are also given 
in this Table. 

k Whether the emission line feature at the red end of [Fe II] A7155, which we identify as [Ni II] A7378 (§4 & 5), is detected or not at least 
at one epoch. 

c The spectra at 358 days and 398 days arc the Mid-IR and NIR ones, respectively. 



and instruments. The signature of the possible asym- 
metry in the deflagration phase, however, can only be 
probed by looking deeper into the innermost regions. In 
this respect, late-phase (~ 1 year past the explosion) 
spectroscopy can provide an important diagnostics. Fol- 
lowing the homologous expansion, the SN ejecta become 
transparent to optical and longer wavelengths, and thus 
emission line profiles can be used to probe the distribu- 
tion of elements that emit the light or that deposit the 
energy. This strategy has been applied to investigate 
the aspherical nature of core-collapse SNe from massive 
stars (e.g., Maeda et al. 2002; Chugai et al. 2005; Maz- 
zali et al. 2005; Maeda et al. 2008; Modjaz et al. 2008; 
Taubenberger et al. 2009; but see also Milisavljevic et 
al. 2009). 

Probing the explosion geometry of SNe Ia using late- 
time spectroscopy is still a fairly young field. Hoflich ct 
al. (2004) presented Near-Infrared (NIR) spectra for SN 
Ia 2003du taken with the Subaru telescope. They dis- 
cussed that a flat-topped profile of the [Fe II] 1.644 /im 
emission feature is likely a result of a hole in the dis- 
tribution of 56 Ni. Motohara et al. (2006) added two 
other Near-Infrared (NIR) spectra for SNe Ia taken with 
the Subaru telescope (SNe 2003hv, and 2005W). The 
[Fe II] 1.257 /im and 1.644 fim emission lines displayed 
different profiles, and especially, those in SN 2003hv were 
blueshifted by - 2,000 - 3,000 km s" 1 . Motohara et 
al. (2006) suggested that this can be interpreted as 
evidence of a global asymmetry in SN 2003hv. Ger- 
ardy et al. (2007) also found a similar blueshift, in the 
[Co III] 11.88 /im line detected in their Mid-Infrared 
(Mid-IR) spectrum of SN 2003hv, taken with the Spitzer 
Space Telescope. On the other hand, the late-phase opti- 
cal spectra, available for more than 20 normal SNe Ia, all 
look very similar to each other, and seem to show little 
evidence for any asymmetry. 

Recently, Leloudas et al. (2009) reported on observa- 



tions of SN 2003hv, including late-time optical spectra. 
Thus, it is now for the first time possible to examine a 
late-time spectrum of a SN Ia all the way from the optical 
to the Mid-IR. We want to investigate if the kinematical 
interpretation for the blueshifts in the NIR and Mid-IR 
forbidden lines of SN 2003hv is consistent with the opti- 
cal spectrum. 

In this paper, we examine late-time spectra of SN 
2003hv, using multi-dimensional radiation transfer cal- 
culations and a simple kinematic model, and arrive at a 
plausible explosion geometry of SN 2003hv. Using the 
structure that was successfully applied to SN 2003hv, we 
then investigate signatures of the asymmetry and the ex- 
pected diversity in the late-phase spectra, resulting from 
various viewing orientations. We also compile published 
late-time optical spectra of SNe Ia, and investigate to 
what extent the expected diversity exists in the data. In 
doing this, we identify a possible signature of the explo- 
sion asymmetry in late-time optical spectra of SNe Ia. 

This paper is organized as follows. In §2, the obser- 
vational data of SN 2003hv and other SNe Ia are sum- 
marized. In §3, our model and the method of multi- 
dimensional spectrum synthesis are presented. The re- 
sults are shown in §4 & 5; §4 focuses on the geometry of 
SN 2003hv, while the implications for SNe Ia in general 
are presented in §5. The paper is closed with a discussion 
and our conclusions in §6. Future observing strategies to 
further constrain the explosion mechanism are also pre- 
sented in §6. 

2. THE SAMPLE OF SN IA LATE- TIME SPECTRA 

The data set compiled for this work are summarized 
in Table 1. For SN 2003hv we focus on the optical spec- 
trum at ~ 320 days after maximum brightness, taken 
with the VLT telescope equipped with FORS1, as re- 
ported by Leloudas et al. (2009). There are two more 
late-time spectra at ~ 110 and 143 days after maximum 
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v [km s 1 ] 

Fig. 1. — The kinematic model explored in the present study. 
The inner ejecta consist of three characteristic zones; the ECAP 
region (red), High-Density (HD) 5B Ni-rich region (blue), and Low- 
Density (LD) 56 Ni-rich region (gray). The model is shown for 
foff.ECAP = 3, 500 km s _1 (the off-set of the ECAP and the HD 
regions) and v a ff = —1, 500 km s~ x (Tab. 2). 

light (Leloudas et al. 2009), which are briefly discussed 
in §6. A NIR spectrum of SN 2003hv was obtained on 
2004 October 6 with the Subaru telescope equipped with 
CISCO/OHS (Motohara et al. 2006), 394 days after 
maximum brightness. Leloudas et al. (2009) generated a 
combined optical-NIR spectrum mapped to 358 days by 
scaling the flux in the spectra mentioned above by inter- 
polating optical and NIR photometric points. We also 
use a Mid-IR spectrum taken by Gerardy et al. (2007), 
with the Spitzer space telescope at 358 days. 

Optical late-time spectra of other SNe la are compiled 
from the S USPECT database 8 . We selected data which 
satisfy two criteria: (1) the spectrum is taken at least 
100 days after maximum brightness, and (2) the spec- 
trum covers the wavelength range ~ 7, 000 - 7, 500 A 
(for this criterion, see §4 & 5). References for the data 
compiled from the SUSPECT database are listed in Ta- 
ble 1. Adding to this, we have also included SN 2000cx 
at the later epoch (+ 360 days) from Sollerman et al. 
(2004), and SN 2001el from Mattila et al. (2005). 

3. MULTI-DIMENSIONAL SPECTRUM SYNTHESIS 
3.1. Model 

Rather than working with multi-dimensional, detailed 
hydrodynamic and nucleosynthesis models, we calculate 
nebular spectra for a simplified kinematic model, and 
try to obtain constraints on the model parameters from 
the late-time observations of SN 2003hv. In the model, 
we have three zones (Fig. 1); the off-set region filled 
with neutron-rich, Fe-peak elements produced with elec- 
tron capture reactions (the ECAP zone), the off-set, rela- 
tively high density region filled with 56 Ni (the HD zone), 
and the extended, relatively low density region filled with 
56 Ni (the LD zone). The model parameters are the fol- 
lowing: 

• icecap: The off-set velocity of the ECAP- and 
the HD-zones, with respect to the zero- velocity cen- 

8 The Online Supernova Spectrum Archive, "SUSPECT", is 
found at http://bruford.nhn.ou.edu/~suspect/indexl.html . 



ter of the ejecta. It is assumed that the off-set ve- 
locity is the same for the ECAP- and the HD-zones 
for simplicity. 

• ^off.LD: The off-set velocity of the LD zone, with 
respect to the zero- velocity center. 

• ^ecapj whd 5 wld= The outer velocities of the 
ECAP-, HD-, and the LD-zones, respectively. The 
distribution of the elements and the density is as- 
sumed to be homogeneous within each region, with 
spherical symmetry with respect to the off-set posi- 
tion (but the HD /LD zone has a hole corresponding 
to the ECAP/HD zone). 

• Mecap» AffiD? Mjjji The mass in each region. 

In each zone, the initial abundance is set as follows; 
In the ECAP-zone, 10% of the mass is in 56 Ni, and the 
remaining 90% is in stable 58 Ni. The other regions are 
assumed to fully consist of 56 Ni. 56 Ni decays to 56 Co 
and then to 56 Fe, and thus the synthetic spectra are a 
mixture of forbidden lines of Ni, Co, and Fe. The masses 
are varied so that the model can roughly reproduce the 
total flux as well as the fluxes of the emission lines of 
interest, as compared with the late-time spectra of SN 
2003hv. The outer velocities are set so that the predicted 
line widths are consistent with the observed widths. The 
off-set velocities affect the wavelength centers of various 
emission lines. 

Note that this model represents only the inner regions 
of the ejecta, and does not include the outer regions (i.e., 
the regions dominated by intermediate mass elements), 
since the outer regions contribute little to the late-time 
emission. Detailed model fit to all the spectral features 
is beyond the scope of this paper, as is evident from the 
simplifications in our model. Also, the ejecta probably 
have a clumpy structure as well in realistic multi-D explo- 
sions, while such an effect is not included in the present 
study (see, e.g., Leloudas et al. 2009). 

The model is constructed to represent the main fea- 
tures of hydrodynamic explosion models, especially of the 
delayed-detonation model. The initial deflagration pro- 
duces neutron-rich Fe-peak elements with electron cap- 
ture reactions, and then 56 Ni as the density decreases at 
the deflagration front (e.g., Nomoto et al. 1984; Iwamoto 
et al. 1999; Brachwitz et al. 2000). We therefore regard 
the ECAP-zone and the HD-zone as the products of the 
early deflagration phase (but see §6 for another scenario 
to create the off-center ECAP-zone). The subsequent 
detonation produces a relatively low-density 56 Ni-rich 
region (the LD-zone). The distribution of the detona- 
tion products is expected to be more or less spherical 
(Ropke & Niemeyer 2007b; Kasen et al. 2009; Maeda et 
al. 2009c). We still allow for a small off-set (opposite 
to the deflagration) for the LD-zone, since the detona- 
tion may well be stronger in the direction where a larger 
amount of fuel is left after the deflagration phase. 

3.2. Multi- Dimensional Spectrum Synthesis 

The input model is mapped onto a three-dimensional 
Cartesian grid, with 51 3 zones. Three-dimensional neb- 
ular spectrum synthesis calculations are then performed, 
using the 7-ray transfer module (Maeda 2006a) and 
the nebular spectrum synthesis module (Maeda et al. 
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TABLE 2 
Model Parameters' 1 



^off,ECAP »off,LD «ECAP 


"HD 


»LD 


AfECAP 




Af L D 


3,500 -1,500 3,000 


5,000 


10,000 


0.1 


0.1 


0.2 



The units for velocities and 



. are km s 1 and Mq , respectively. 
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Fig. 2. — Synthetic late-time spectra of the SN la models, as 
compared to the observed spectrum at day +320 (gray: dereddened 
with E(B — V) = 0.016 mag and Ry = 3.1, the redshift corrected 
for the host galaxy recession velocity, and the flux calibrated to 
day +358. The model spectrum is calculated for 375 days after 
the explosion, and the flux is converted to the observed flux using 
the distance modulus fi = 31.37. The W7 model of Nomoto et al. 
(1984) is a spherically symmetric model, and does not include any 
kinematic off-set. The model fluxes are shown in logarithmic scale, 
with offsets of -2 (W7: black), -4 (off-center, deg: red), and -6 
(off-center, 180 deg: blue). 

2006c) of the SAMURAI code - the SupernovA Multi- 
dimensional Radiation transfer code 9 . As a reference, 
we have also performed one-dimensional nebular spec- 
trum synthesis for the classical, deflagration model W7 
(Nomoto et al. 1984). 

Details of the calculation method are given in Maeda 
(2006a) (for 7-ray transport), Mazzali et al. (2001) 
(for one-zone nebular spectrum synthesis), and Maeda 
et al. (2006c) (for multi-dimensional spectrum synthe- 
sis). Good reviews on the nebular spectrum synthe- 
sis are given by e.g., Axelrod (1980), Ruiz-Lapuente & 
Lucy (1992), Kozma & Fransson (1998ab), and Liu et al. 
(1998). 

The 7-rays emitted by the decay chain 56 Ni — + 56 Co — > 
56 Fe produce non-thermal high energy electrons mainly 
by Compton scatterings, and the electrons give rise to 
impact ionization and excitation. The energy is also lost 
in thermalization and heating of the ejecta, through in- 
teractions with thermal electrons. An additional energy 
input comes from positrons emitted by the radioactive 
decays, with a stopping length much shorter than 7-rays. 
In our calculations, it is assumed that the positrons are 
fully trapped on the spot. The deposition of 7-rays and 

9 The SAMURAI is a compilation of 3D codes adopting Monte- 
Carlo methods to compute the high-energy light curve and spec- 
tra (Maeda 2006a), optical bolometric light curve (Maeda et al. 
2006b), and optical spectra from early (Tanaka et al. 2006, 2007) 
to late phases (Maeda et al. 2006c). 



Fig. 3. — Same as Figure 2, but for NIR wavelengths. The 
original observed spectrum (Motohara et al. 2006) is again flux- 
scaled to day +358, and compared to the synthetic spectra at 375 
days after the explosion. The model fluxes are shifted by the same 
amounts as in Figure 2. 

non-thermal positrons is the dominant heating process 
in the SN ejecta. 

At late phases, the heating is balanced by cooling 
from various thermal-collisionally excited emission lines, 
mostly forbidden. The non-thermal excitation is ne- 
glected in our calculations. The non-thermal impact 
ionization is balanced with recombination. As the den- 
sity is low following the expansion, and the energy in- 
put is not dominated by thermal photons, the result- 
ing ionization stages and level populations are gener- 
ally not in Local Thermodynamic Equilibrium (LTE). 
For a given (thermal) electron temperature (T ), the 
ionization-recombination balance is solved (without pho- 
toionization in our calculations; see Sollerman et al. 
2004), providing the electron density (N e ) and ionization 
stages for each element. Once the ionization stage and 
the electron density are given, detailed balances for level 
populations are solved for each ion, under the condition 
of global heating-cooling balance. The balance then pro- 
vides the electron temperature, as well as the resulting 
emission-line spectrum. In our numerical calculations, 
the ionization balance and the detailed balance for level 
populations are solved iteratively, until the electron den- 
sity and temperature converge at every mesh point. 

The spectrum is sampled into 19 zenith angles, divided 
by 10 degrees each. The wavelength resolution in sam- 
pling the final spectrum corresponds to ~ 450 km s _1 , 
roughly the same width as the spatial resolution of the 
input model. Note that any structure in the synthetic 
spectrum below this resolution is not reliable. 

4. SN 2003HV 
4.1. Model Overview 
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TABLE 3 
Typical Values in Model Output 3 



-0.03 0.00 0.03 

v/c 



-0.03 0.00 0.03 



Fig. 4. — Comparison between synthetic and observed line pro- 
files. The observations are shown in gray (see Table 1 for refer- 
ences). The black line is the W7 model artificially shifted to the 
blue by 2,800 km s —1 . Our off-set model is shown for the observers 
in the direction of the off-set of the ECAP-zone (red) and in the 
opposite direction (blue). Unlike the models in Figures 2 and 3, 
the observed profiles are compared to the synthetic spectra at 333 
days after the explosion in the optical (day +320) and at 410 days 
after the explosion in NIR (day + 398), in order to avoid possi- 
ble evolutionary effects in the line profiles. The flux is arbitrarily 
normalized in this comparison. 



Figures 2 and 3 show synthetic spectra for the W7 
and our off-set model (Fig. 1), in the optical (Fig. 2) 
and in the NIR (Fig. 3), at 375 days after the explosion. 
In these figures, we show the observed spectrum at 320 
days (optical) and 398 days (NIR) after the maximum, 
but with the flux calibrated to 358 days by Leloudas et 
al. (2009). Figure 4 shows the comparison of line profiles 
for some selected emission lines. The model parameters 
are summarized in Table 2. 

Despite the simplification, the kinematic model can ac- 
count for the presence of the major emission lines in the 
late-phase spectrum of SN 2003hv, and therefore also for 
those of SNe la in general, as do the W7 model (see 
Axelrod 1980, who showed that the main features of a 
SN la nebular spectrum can be modeled by a simple, 
homogeneous 56 Ni-rich sphere). The observed flux, es- 
pecially in the NIR, is smaller than that predicted by 
the W7 model, and consistent with that in the kinematic 
model. This means that the mass of 56 Ni is ~ O.3M0 
in SN 2003hv, smaller than the typical value in SNe la 
(~ 0.5 - 0.6M Q , e.g., Mazzali et al. 2007) and than the 
W7 value (~ 0.6M Q ). This value is roughly consistent 
with that estimated from the peak luminosity (~ 0.4M Q ) 
(Leloudas et al. 2009). We note that this is also re- 



Zone 


T c 


logiVo 


Fe+ Fraction 13 


ECAP 


2,000 


5.7 


90% 


HD 


7,000 


5.4 


30% 


LD 


11,000 


4.8 


5% 



a Obtaincd as a result of spectrum synthesis 
for 375 days after the explosion. The units 
for T e (the electron temperature) and N G (the 
electron density) arc in cgs. 

^The fraction of Fe^. The remaining frac- 
tion is in Fe + + . 

lated to relative contributions from different zones in our 
model, indicating that the W7 model has too much high 
density, low temperature materials to be consistent with 
the NIR flux. 

In Figures 2 and 3, the expected rest wavelengths of 
some selected emission lines are marked by lines: [Fe III] 
blend at 4700A, [Fe II] A7155, and [Ni II] A7378 in the 
optical (Fig. 2), [Fe II] 1.257/um and [Fe II] 1.644/xm in 
the NIR (Fig. 3). As found by Motohara et al. (2006), 
the observed NIR emission features are blueshifted as 
compared to those predicted by a spherically symmet- 
ric model like W7. Leloudas et al. (2009) pointed out 
that [Fe II] A7155 and [Fe II] A8621 also show a similar 
amount of blueshift as compared to the expected wave- 
length. Furthermore, they suggested that [Fe II] A8621 
can be one of the best and cleanest lines to follow the 
geometry of the explosion. 

We here suggest that the feature at ~ 7300 A is dom- 
inated by [Ni II] A7378, and then note that this feature 
is also blueshifted as compared to the expected position. 
At the same time, features in the blue part of the optical 
spectrum, e.g., at ~ 4700 A and ~ 5250 A, do not show 
the blueshifts. 

This is more clearly seen in Figure 4, where the syn- 
thetic spectrum for the W7 model is artificially shifted 
in wavelength, by 2, 800 km s -1 to the blue, to fit the 
NIR features (see Motohara et al. 2006). The blueshift 
seen in [Fe II] A7155 is also explained by the same model. 
However, the features in the blue (~ 4700, 5250 A) should 
also shift to the blue, and no longer fit the observed wave- 
lengths. 

The inconsistency in the wavelengths is a strong ar- 
gument against a simple, bulk off-set model as the 
origin of the wavelength shift in the NIR (and at 
7155,7378,8621 A). Our 2D model (though computed 
in 3D) was constructed to overcome this inconsistency. 
As shown in Figures 2-4 (especially evident in Figure 
4), the observed wavelengths of strong lines, in the opti- 
cal, NIR, and Mid-IR wavelengths are all well explained 
by our model, with the viewing orientation close to the 
direction of the off-set. 

4.2. Ionization /Thermal Structure and Line Shift 

These model results can be understood on the basis 
of the ionization and thermal structure. Table 3 shows 
the typical electron temperature, electron density, and 
ionization stage in each of the zones. The ECAP-zone is 
in a low ionization state (mostly singly-ionized species), 
because of the small energy input (inefficient ionization) 
and high density (efficient recombination). Despite the 
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TABLE 4 
Optical and Infra-red Lines 3, 



XA/tlWoloVl ITT f\ ill YY°i i 

VVcLVL-lCIl^LIl yfJjlll J 


Ion 
on 


Term 


F (rTn _1 ) b 


Shift c 




0.4658 


Fe III 


5 D 4 - |F 4 
5 D 3 - 3 F 3 


21462.2 


No 


LD 


0.4701 


Fe III 


21699.9 


No 


LD 


0.4734 


Fe III 


5 D 2 - |F 2 


21857.2 


No 


LD 


0.5262 


Fe II 


a 4 F 7/2 - a 4 H 11/2 


21430.4 


No 


LD 


0.7155 


Fe II 


a 4 Fg/ 2 - a 2 Gg/2 


15844.7 


Yes 


HD 


0.7378 


Ni II 


2 D 5/2 - 2 F 7/2 


13550.4 


Yes 


ECAP 


0.8617 


Fe II 


a 4 F 9/2 - a 4 P 5 / 2 


13474.4 


Yes 


HD 


1.257 


Fe II 


a 6 Dg/ 2 - a 4 D 7 / 2 


7955.3 


Yes 


HD 


1.644 


Fe II 


a 4 Fg/ 2 - a 4 D 7/2 


7955.3 


Yes 


HD 


2.218 


Fe III 


3 H 6 - 3 G 5 


24558.8 


No 


LD 


2.348 


Fe III 


3 H 5 - 3 G 5 


24558.8 


No 


LD 


2.874 


Fe III 


?F 3 - 3 G 3 
|F 3 - 3 G 3 


24940.9 


No 


LD 


2.904 


Fe III 


25142.4 


No 


LD 


3.228 


Fe III 


|f 4 - 3 g 5 


24558.8 


No 


LD 


4.114 


Fe II 


a 6 D 9/2 - a 4 F 7/2 


2430.1 


Yes 


HD + ECAP 


4.606 


Fe II 


a 6 D 5/2 - a 4 F 5/2 


2838.0 


Yes 


HD + ECAP 


4.888 


Fe II 


a 6 D 1/2 a 4 F 3/2 


2430.1 


Yes 


HD + ECAP 


5.339 


Fc II 


a 6 D 9/2 - a 4 F 9/2 


1872.6 


Yes 


HD + ECAP 


6.634 


Ni II 


2 D 5/2 - 2 D 3/2 


1506.9 


Yes 


ECAP 


7.350 


Ni III 


3 F 4 - 3 F 3 


1360.7 


Yes 


ECAP 


10.52 


Co II 


a 3 F 4 -a 3 F 3 


950.51 


Yes 


HD 


10.67 


Ni II 


4 Fg/ 2 - 4 F 7/ / 2 


9330.0 


Yes 


ECAP 


11.00 


Ni III 


3 F 3 - 3 F 2 


2269.6 


Yes 


ECAP 


11.88 


Co III 


a 4 F 9/2 - a 4 F 7/2 


841.2 


Yes 


HD + LD 


12.72 


Ni II 


4 F 7/2 - 4 F 5/2 


10115.7 


Yes 


ECAP 


14.74 


Co II 


a 5 F 5 -a 5 F 4 


4029.0 


Yes 


HD 


17.93 


Fc II 


a 4 F 9/2 - a 4 F 7/2 


2430.1 


Yes 


HD + ECAP 



a This is not a complete line list; especially, in the optical here are 
main text. The NIR and Mid-IR line list covers the strongest lines 
^The upper energy level of the transition. 
c Shift in the wavelength depending on the viewing angle. 
^Region which makes the predominant contribution. 



only the lines discussed in the 
in the model. 



low ionization state, the electron density is high (because 
of large material density). The electron temperature is 
accordingly low. In the LD-zone, the situation is the 
opposite. It is in a high ionization state (mostly doubly- 
ionized), the electron density is small and thus the elec- 
tron temperature is kept as high as ~ 10, 000 K. The 
HD-zone has properties intermediate between these two 
zones. 

The characteristic strong emission lines in the syn- 
thetic spectrum are summarized in Table 4. The strong 
feature at <~ 4, 700A is a blend of several [Fe III] lines, all 
of which are formed under the high ionization and high 
temperature condition. The feature is thus dominated 
by emissions from the extended LD zone, which is dis- 
tributed in an approximately spherical way. Therefore, 
the feature at <~ 4700A does not show a large velocity 
shift. The [Fe II] 1.257 fiia and 1.644 fim, on the other 
hand, come from the low ionization and low tempera- 
ture conditions, and thus are emitted mostly from the 
off-set HD-zone. Therefore, these NIR [Fe II] lines show 
a large blueshift if viewed from the direction of the off- 
set. [Ni II] A7378 is emitted from the ECAP zone (i.e., 
58 Ni), and show a large velocity shift. Contribution from 



radioacti 



3 Ni to this line is negligible at these late 



phases because of the short decay time scale. 

The behavior of other [Fe II] lines in the optical range 
can be understood in the same manner, but is compli- 
cated by a competition between the ionization and tem- 
perature effects. [Fe II] AA7155 and 8621 show the veloc- 



ity shift, as these are mostly emitted from the HD-zone. 
This is due to the low ionization in this zone. On the 
other hand, the feature at ~ 5250 A does not show a 
large shift (Fig. 2), which is not interpretable by the ion- 
ization effect alone. The strongest line in this feature is 
[Fe II] A5262. The excitation temperature for this line is 
high, and thus the high temperature in the LD-zone is 
preferred, despite the small amount of Fe + present there. 
Also, the contribution of [Fe III] A5270 from the LD-zone 
is not negligible. As a result, the large contribution to 
the emission feature at <~ 5250 A comes from the LD- 
zone. This explains why this feature does not show a 
large velocity shift. 

[Co III] 11.88 /xm is a ground state transition with an 
excitation temperature of only ~ 1,000 K, much lower 
than the temperatures in either the HD or LD-zone. In 
addition, the fraction of the "doubly" -ionized Co is sim- 
ilar between these two zones (unlike the singly-ionized 
Co). The LD and HD zones thus provide comparable 
contributions to the [Co III] 11.88 feature. As a re- 
sult, the line profile is a combination of a broad compo- 
nent centered at the rest wavelength and a relatively nar- 
row component whose central wavelength is blucshiftcd 
if viewed from the off-set direction. Because of the strong 
narrow component, the line as a whole shows a velocity 
shift, depending on the viewing direction. 



SIGNATURES OF EXPLOSION ASYMMETRY IN SNE IA 
IN GENERAL 
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Fig. 5.— Analysis of the [Ni II] A7378 line profiles in the 12 
SNe la. The velocity is set assuming that the rest wavelength is 
at 7378 A. (a) Observed line profiles. The rest wavelengths of 
[Fe II] A7155 and [Ni II] A7378 are shown by dotted lines, (b) 
[Ni II] A7378 in observations, after removing the underlying con- 
tinuum (or possible other lines) as much as possible (see main text), 
(c) Synthetic line profiles of the [Ni II], depending on the viewing 
orientation. 

In §4, we have shown that the somewhat puzzling fea- 
tures in the late-time spectra of SN 2003hv, i.e., some 
lines showing a blueshift while others showing no veloc- 
ity shift, can be explained by a geometrical effect. In this 
section, we expand the discussion in §4; we investigate 
the explosion geometry of SNe la in general. We also 
present diagnostics of the geometry using NIR and Mid- 
IR lines (Tab. 4), for existing as well as for future obser- 
vatories such as JWST (James Webb Space Telescope) 
or SPICA (SPace Infrared telescope for Cosmology and 
Astrophysics mission). 

5.1. [Ni II] X7378 as diagnostics of the deflagration 

phase 

It has been believed that signatures of ejecta asymme- 
try are not evident in the optical range (§1 and §5.2). 
Figure 5 shows the feature around ~ 7000 — 7500 A, i.e., 
[Fe II] A7155 (with some contribution from [Fe II] A7171) 
and [Ni II] A7378, for the 12 SNe la (Table 1; see below). 

Contrary to the earlier expectation, we find that 
there is indeed a probable signature of ejecta asymme- 
try. As shown in Figure 5, the feature shows a varia- 
tion in the central wavelength (after removing the host 
galaxy redshift). In many cases, the shift is larger than 
1,000 km s _1 , which is too much to be due to a pecu- 
liar motion of the SN progenitor with respect to the host 
center. The [Fe II] A7155 and [Ni II] A7378 always show 
a similar degree of the wavelength shift. 

Both blucshifts and redshifts exist, and the shift does 
not appear to correlate with the age of the SN (see §6 
for further discussion). These observed characteristics 
suggest that the shift is not caused either by radiation 
transfer effects or by other unidentified emission lines. 
Thus, the mounting evidence is that the variation of the 
central wavelengths represents a real variation in the line- 
of-sight velocity of the region emitting [Fe II] A7155 and 
[Ni II] A7378. 

As shown in §4, these lines are emitted from deep 
parts of the ejecta: [Fe II] A7155 from the HD-zone and 
[Ni II] A7378 from the ECAP-zone in our model. We 
point out that especially the [Ni II] should preserve im- 
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Fig. 6.— Analysis of the [Fe III] blend at 4700A in the 12 SNe 
la. The velocity is set assuming the rest wavelength of 4700 A. (a) 
Observed line profiles. The rest wavelengths of [Fe III] A4701 is 
shown by a line, (b) Synthetic line profiles of the [Fe III] feature, 
depending on the viewing orientation. 

portant information on the explosion dynamics, since the 
ECAP-zone is attributed to the region created by the ini- 
tial deflagration (or, the detonation passing through the 
central region of the white dwarf if the initial deflagra- 
tion is very weak; §6), with efficient electron capture re- 
actions. This is supported by the relatively narrow width 
(< 3, 000 km s" 1 ) of the [Ni II] line. 

One may argue that the central wavelength of the 
[Ni II] A7378 could apparently shift to the blue because 
the blue wing of this feature is contaminated by the 
[Fe II] (and by [Ca II] AA7291, 7307, although it does 
not likely contribute much 10 ). To check this, we tried 
to obtain intrinsic profiles of the [Ni II] line by subtract- 
ing the possible contamination of other lines, especially 
in the blue wing (Fig. 5). Practically, we assumed a 
(pseudo)-continuum connecting the flux minima on both 
side of 7, 380 A (the blue minimum corresponds to the 
valley between the [Fe II] and [Ni II]). This continuum 
flux generally decreases towards the red, and thus sub- 
tracting it could result in a shift of the peak wavelength 
to the red to some extent. However, we found that this 
effect is small, as compared to the observed variation. 

Another possible concern is the simplified abundance 
distribution in our model, in which the HD and LD-zones 
are assumed to contain no stable Ni. Even without elec- 
tron captures realized in the early deflagration phase, 
some amount of 58 Ni should be produced in the region 
where 56 Ni is the dominant burning product. The typ- 
ical mass fraction of 58 Ni in such a region is ~ 5% for 
solar metallicity (e.g., Iwamoto et al. 1999; Timmes et 
al. 2003). To check the effect, we replaced 5% of the 
material by 58 Ni in the model and repeated the calcula- 
tions. We thereby confirmed that the result is not de- 
feated by this change, and that [Ni II] A7378 emitted 
from the HD /LD-zones is negligible as compared to that 
from the ECAP-zone. 

The [Ni II] A7378 line is not always strong in the ob- 
served spectra of SNe la. We find that 12 SNe la (Fig. 5) 

10 We have not included Ca in our model. The W7 model pre- 
dicts that the contribution from [Ca II] A7291, 7307 is about 10% 
(Leloudas et al. 2009). Also, the W7 model predicts broad [Ca II] 
which is not compatible to the observed relatively narrow emission 
features. 
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Fig. 7.— Analysis of NIR [Fe II] line profiles in the 3 SNe la 
(Roflich et al. 2004; Motohara et al. 2006). The velocity is set 
assuming the rest wavelength of 1.257 fira (J-band) and 1.6440 (im 
(H-band). (a) Observed line profiles, (b) Synthetic line profiles of 
[Fe II] 1.257 Atm. (c) Synthetic line profiles of [Fe II] 1.644 /tm. 

show this feature among the 20 SNe we investigated (Ta- 
ble 1). The number of SNe la, showing the [Ni II] A7378 
stronger than the [Fe II] A7155 is even smaller. Although 
[Ni II] A7378 is mainly emitted from the ECAP-zone 
in our present calculations, the contribution from the 
HD/LD regions may not be negligible in case of larger 
Mhd and/or Mld- Also, with an increasing amount of 
56 Ni, [Fe II] A7388 and [Fe II] A7452, could hide the [Ni II] 
feature. Thus, the observed diversity in the detection of 
the [Ni II] may indicate that Mecap / (A/hd+^ld) varies 
among objects (e.g., Mazzali et al. 2007). 

5.2. [Fe III] blend as diagnostics of the outer region 

In Figure 6, we show line profiles of the same sample of 
SNe la, but centered at 4700 A. Unlike the [Ni II] A7378, 
there is no large variation. The absence of a significant 
wavelength shift has been a strong argument against any 
global asymmetry in SNe la, since this feature is the 
strongest in SN la nebular spectra and has thus naturally 
been used to infer the geometry. 

This feature is a blend of several [Fe III] lines, with 
the wavelength separations smaller than the typical line 
width (10,000 km s" 1 ); [Fe III] A4658, [Fe III] A4701, 
[Fe III] A4734, [Fe III] A4755, [Fe III] A4769, and 
[Fe III] A4778. As a result of the blending, the central 
wavelength of the feature is ~ 4, 700A. Note that this 
central wavelength is not sensitive to underlying models 
as long as spherically symmetry is assumed, since the 
excitation temperature of these lines are all similar and 
thus the relative contribution is basically determined by 
the transition probabilities. 

With the central wavelength of 4, 700A, no redshifts are 
observed in these SNe la. There are, on the other hand, 
some SNe la showing small blueshifts in this feature. 
This argues against the idea that the wavelength shifts 
here are caused by a geometrical effect. We note that the 
SNe la that show relatively large blueshifts (> 1, 500 km 
s _1 ) are mostly young objects (the spectra taken before 
day +200) - SNe 2004dt (+152), 1990N (+186), 1994D 
(+106), 2000cx (+147), and 2001V (+106). The cen- 
tral wavelength of the 4,700A feature in spectra taken 
at > 200 days is clustered at the expected wavelength 

(~ 4,700 A). This is further discussed in §6. We con- 
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Fig. 8. — Analysis of Mid-IR line profiles. Synthetic line profiles 
are shown for various lines (a — f; see Tab. 4). Three lines in 
each panel are for different viewing angles (0, 90, 180° from top to 
bottom; red, green, and blue, respectively). 

elude that the wavelength shift in the 4700 A feature 
is not caused by geometric effects, but by either a radi- 
ation transfer effect or contamination from other lines, 
e.g., Mg I] A4571. 

The behavior that the [Fe III] blend is always at ~ 
4,700A in sufficiently late phases (> +200 days) can be 
understood from our model results. Since this blend is 
emitted from the LD zone, it indicates that the LD zone 
is more spherically distributed than the inner ECAP and 
HD zones. 

5.3. NIR [Fe II] lines as a promising tool to investigate 
the geometry 

Figure 7 shows a similar analysis for the NIR (J and 
H) [Fe II] features - [Fe II] 1.257 /im and 1.644 /mi. As 
discussed in §4, these lines are mostly emitted from the 
HD-zone, and thus show a large variation as a function 
of the viewing angles. 

Not many observational data are available for these 
NIR lines, because of the difficulty raised by the OH 
airglow lines. In Figure 7 we show the data presented by 
Motohara et al. (2006). NIR spectra for two other SNe la 
have been reported: SN 1991T (Bowers et al. 1997) and 
SN 1998bu (Spyromilio et al. 2004). The [Fe II] 1.644 /im 
line in SN 1991T is almost symmetric with respect to the 
rest frame of the host galaxy, while SN 1998bu seems to 
show a small degree of blueshift. 

In principle, it is favorable to check the consistency 
of the wavelengths of various lines, as we have done for 
SN 2003hv. [Fe II] 1.644 /im in SN 2003du seems to 
show some blueshift (Hoflich et al. 2004; Motohara et 
al. 2006), although the Signal-to-Noise ratio (S/N) is 
not good. [Ni II] A7378 in SN 2003du is only marginally 
detected (Fig. 5), and the poor S/N makes it difficult to 
determine the exact position of the center of the line. In 
SN 1998bu, [Ni II] A7378 is relatively well identified with 
better S/N than in SN 2003du, and it shows - 1, 000 km 
s _1 of the blueshift, at least qualitatively in agreement 
with the NIR observation. 

5.4. Mid-IR lines as a new probe 

In the mid-IR wavelength range, there are various iso- 
lated forbidden lines of Fe, Co, and Ni. There are also 
lines from intermediate mass elements such as [Ar II] and 
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Fig. 9. — (a) The central wavelength of [Ni II] A7378 (red filled squares) and the [Fe III] blend at 4700 A (black filled squares) in the 
synthetic spectra as a function of the viewing angle. Note that the binning in the synthetic spectrum is about 0.0015c, and fluctuations at 
this level can simply be numerical noise. The angle dependence of [Ni II] A7378 is consistent with v/c ~ O.O12cos(0) as expected. On the 
other hand, the 4700A feature is not sensitive to the viewing angle. Panels (b) and (c) show the line shift derived for the 12 SNe la. The 
SNe are divided according to their phase: (b) SNe after day +200, and (c) SNe before day +200. The red symbols are for [Ni II] A7378, 
while the black ones are for the [Fe III] blend at 4700A. For both model and observed spectra, the central wavelength is defined as the 
wavelength with 50 % of the total line luminosity on either side. 



[Ar III] (Gerardy et al. 2007), which are not included in 
our model. For these longer wavelengths, lines tend to be 
more isolated, and thus hardly affected by a line blend- 
ing. The lines in the mid-IR therefore provide a direct 
view on the geometry of the emitting region. 

In Figure 8, some selected lines synthesized in our 
model are shown. Emission features at ~ 2.2 /an and 
~ 2.9 /im are both dominated by [Fe III] with an ex- 
citation temperature of ~ 30, 000 K. Thus, these lines 
behave in the same way as the [Fe III] blend at 4,700A; 
virtually no wavelength shift is seen, irrespective of the 
viewing orientation. 

The behavior of [Ni II] 6.634 /im is basically the same 
as in [Ni II] A7378. This line provides a very strong 
diagnostics on the geometry of the initial deflagration 
phase, even better than [Ni II] A7378: (1) the excitation 
temperature is only a few 1,000 K (thus no doubt it is 
emitted from the ECAP-zone where not much heating 
is taking place), and (2) there is no blending with other 
lines. These effects are already taken into account in the 
analysis of [Ni II] A7378, but [Ni II] 6.634 /im can provide 
the the same diagnostics in a more model-independent 
way. 

On the other hand, the analysis of [Ni III] 7.350 /jm 
may be complicated by the competition between the high 
ionization and low excitation temperature. In our model, 
it is dominated by the emission from the ECAP-zone. We 
have confirmed that adding a representative amount of 
58 Ni in the HD/LD-zones does not affect our result (§5). 

[Co III] 11.88 fim has already been discussed in §4, 
in the context of SN 2003hv. From Figure 8, it is seen 
that this line has two components; a broad symmetric 
feature arising from the LD zone and a narrow com- 



ponent from the HD-zone. The former is always at 
the rest wavelength, but the latter shows a variation in 
the central wavelength as a function of the viewing an- 
gle. The relative strength may well vary depending on 
Mhd/Mld- In the present model, the ratio is set by 
the requirement that the fluxes in the blue part of the 
optical range (mainly emitted from the LD-zone) and in 
the NIR (mostly from the HD zone) should be consistent 
with the observed spectra of SN 2003hv. 

[Fe II] 17.93 /im behaves in a way similar as 
[Fe II] 1.644 fim. A large variation in the central wave- 
length is seen. One important difference is that we do see 
a contribution from the ECAP-zone. Because of the low 
excitation temperature of [Fe II] 17.93 /im, this comes 
from a small amount of the decay product of 56 Ni in the 
ECAP-zone where the electron temperature is low. As 
a result, this line does not show a flat-topped profile, as 
seen in [Fe II] 1.644 /im (§6 for further discussion). 

6. DISCUSSION AND CONCLUSIONS 

6.1. Summary 

In this paper, we have investigated the geometry of 
the innermost region of SNe la, which should provide 
information on the explosion mechanism. We have shown 
that some enigmatic properties observed in the late-time 
optical, NIR, and Mid-IR spectra of SN la 2003hv - that 
some emission lines show blueshift while others do not - 
can in fact be naturally explained by a simple kinematic 
model, if the viewing angle is close to the direction of 
the off-set. In this model, the high density regions (the 
ECAP-zone dominated by stable Fe-peak elements, and 
the HD-zone dominated by 56 Ni) are off-set with respect 
to the center of the expansion by ~ 3,500 km s _1 , and 



10 



0.4 



'■5 0.2 



Obs. 


§ 


/ 


Expectation 




































































n 


11 





































-0.010-0.005 0.000 0.005 0.010 
v/c 



Fig. 10. — The number distribution of SNe la as a function of 
the central wavelength of the [Ni II] A7378 feature. The observed 
number distribution, using 12 SNe la, is shown by shaded bars. 
The expected distribution from the present model (as a function 
of the viewing angle) is shown by the line. 

are surrounded by a less asymmetric, low density 56 Ni- 
rich region (the LD-zone). 

Table 4 summarizes the expected qualitative behavior 
of strong emission lines in late-time nebular SN la spec- 
tra. Forbidden lines from ions with high ionization stage 
and/or high excitation temperature are mostly emitted 
from the outer, "spherical" LD-zone. Such lines show lit- 
tle variation in the central wavelength, irrespective of the 
viewing angle. On the other hand, lines from ions with 
low ionization stage and/or low excitation temperature 
are generally emitted from the inner, "off-set" ECAP- 
or HD-zones. These lines show a strong variation in the 
central wavelength as a function of the viewing angle. 

Encouraged by this finding, we have investigated avail- 
able optical late-time spectra of SNe la compiled mostly 
from the SUSPECT archive. Figure 9 summarizes our re- 
sults from the model calculations, and the measurements 
in the 12 observed SNe la sample. From the model, 
we see that the variation in the central wavelength of 
[Ni II] A7378 basically follows the simple kinematic ex- 
pectation. Note, however, for a precise prediction of the 
behaviors of emission lines, relative contributions from 
the different zones have to be considered, which further 
depend on the ionization structure, the electron temper- 
ature and the density. For [Ni II] A7378, the effects of 
these details turn out to be relatively small, but some 
lines are indeed affected by these complications (§5). We 
have provided an explanation on the observational be- 
havior that [Ni II] A7378 shows the diversity in the cen- 
tral wavelength while [Fe III] blend at 4,700A does not: 
they trace different zones, and the inner ECAP and HD 
zones have the large off-set while the outer LD zone does 
not. 

We suggest that [Ni II] A7378 to be one of the few 
lines at optical wavelengths to provide strong diagnos- 
tics of the innermost region of SNe la. Among the 20 
SNe la, about half of them (12 SNe) show a probable 
detection of this feature. Similar diagnostics can be ob- 
tained by [Fe II] A7155 as well. Thus, the spectra in 
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Fig. 11. — Examples of the temporal evolution in (a) the re- 
gion around the [Fe III] blend at 4700A and (b) that around 
[Ni II] A7378. The data are taken from Cristiani et al. (1992) 
for SN 1986G, Leloudas et al. (2009) for SN 2003hv, and Li et al. 
(2001) and Sollerman et al. (2004) for SN 2000cx. 

the wavelength range of - 7, 000 - 7, 500A should pro- 
vide useful diagnostics of the explosion asymmetry. With 
present instrumentation, these features can be more eas- 
ily observed than other strong indicators in the NIR and 
Mid-IR regime. A late-time spectrum of an overluminous 
SN la 2006gz (Hicken et al. 2007) also shows a hint of the 
possible blueshift in these features (Maeda et al. 2009b). 
However, the spectrum is very noisy and the identifica- 
tion is not secure. Also, the spectrum was peculiar in 
a sense that it did not show the strong features in the 
blue (i.e., at ~ 4, 700A). For these reasons, we have not 
included SN 2006gz in our analysis, despite the possible 
importance of asphericity in this object (Hillebrandt et 
al. 2007; Sim et al. 2007; Maeda & Iwamoto 2009a). 

The geometry we derived could naturally arise in SN 
la explosions. The initial deflagration may proceed in a 
very asymmetric way, while the subsequent detonation is 
expected to leave 56 Ni more spherically distributed. An 
important question is whether the geometry suggested 
for SN 2003hv represents a special case or a general struc- 
ture of SN la explosions. Figure 10 shows the distribu- 
tion of 12 SNe la, as a function of the velocity shift in 
[Ni II] A7378. Although the sample is small, the general 
agreement between the observations and the expectation 
assuming the geometry of SN 2003hv as a generic feature 
of SN la explosions, is very encouraging. From this anal- 
ysis, it seems that SN 2003hv may not be very different 
from other SNe la, and its largest velocity shifts can be 
attributed to the viewing angle effect. 

6.2. Temporal Evolution 

We have found that the line center of [Ni II] A7378 does 
not correlate with the epoch at which the spectrum was 
taken. Therefore, the line velocity shift should reflect the 
intrinsic geometry and variation in the viewing angle. On 
the contrary, the central wavelength of the [Fe III] blend 
at 4700A does show a clear correlation with the epoch; 
The line always shows a blueshift before ~ 200 days, 
while there is no significant wavelength shift in spectra 
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taken after ~ 200 days (Fig. 9). This is exemplified by 
the temporal evolution of the nebular spectra, available 
for a few SNe la (Fig. 11): Thus, we suggest that the 
[Fe III] traces the intrinsic geometry only after ~ 200 
days, and that the outer, relatively low density region is 
less aspherical than the innermost region. 

Figure 11 shows that there is also diversity in the tem- 
poral evolution of [Ni II] A7378. In SNe 1986G and 
2003hv, [Ni II] A7378 persisted over a long period, while 
it only shows a transient feature in SN 2000cx. This may 
indeed be consistent with the model calculation; in our 
calculation, the strength of [Ni II] A7378 is largely reduced 
at ~ 350 days (Figs. 2 & 4), following the temperature 
drop and the so-called infrared catastrophe in the ECAP 
zone (see also §6.3). 

The central wavelength of the [Ni II] line does not 
evolve in SN 1986G, while it shows a somewhat larger 
blueshift in the later phase in SN 2003hv. The difference 
is likely related to the distribution of different species, 
and the resulting 7-ray deposition and potentially the 
positron escape. An increasing degree of the blueshift 
in SN 2003hv as a function of time might indicate that 
56 Ni is distributed inhomogeneously within the ECAP 
zone, being relatively concentrated at the high velocity 
"edge" of this zone (e.g., - 2,000 - 3,500 km s" 1 ). In 
such a scenario, 7-rays irradiate the whole ECAP zone 
early on, while later the positrons heat only the relatively 
" high- velocity" region. 

The number of SNe la with a good temporal coverage 
at late phases is still small, and it is strongly encouraged 
to perform multi-epoch spectroscopy. Such intensive ob- 
servations are useful to understand (1) the detailed distri- 
bution of 56 Ni, and (2) detailed thermal processes (e.g., 
7-ray deposition, positron escape, IR catastrophe) in the 
SN la nebulae. It is also important to investigate (3) 
a possible evolution in the line profiles if we are to use 
them to investigate the details of the explosion physics 
(see §6.3 for more details). 

6.3. Future Perspectives 

The conclusion in the present paper that the late-time 
nebular spectra can be used to investigate the geome- 
try of ejecta, thus the explosion mechanism of SNe la, 
provides a new strategy to clarify the nature of SN la 
explosions. Presently, the sample is still limited: The 
late-time nebular spectra are available for ~ 20 SNe la 
in the optical wavelength. There are only a few examples 
in the NIR and Mid-IR. For most of them, a temporal 
sequence of late-time spectra are not available. In this 
section, we highlight the importance of expanding the 
sample of late-time nebular spectra of SNe la. 

In doing this, we summarize predictions from recent 
theoretical models. Two scenarios have been proposed 
to create the ECAP zone in the context of the asymmet- 
ric ignition of the initial deflagration bubbles. In the first 
scenario, many ignition points are clustered near the cen- 
ter of the white-dwarf with possible bulk off-set (§1), pro- 
ducing the ECAP zone by the deflagration flame. In the 
second scenario, one or a few ignition points are placed 
far away from the center of the white dwarf. The defla- 
gration is weak and does not produce much neutron-rich 
materials. The deflagration flame rises to the surface 
without initiating the detonation. The detonation may 
then be triggered near the white dwarf surface, at the op- 



posite side of the initial deflagration ignition (so-called 
Gravitationally Confined Detonation model; Jordan et 
al. 2008). In this case, the initial deflagration is so weak 
and the white dwarf suffers from virtually no expansion 
before the detonation passes through the central region. 
Thus, the electron capture reactions in the detonation 
are not negligible, producing the ECAP zone (Meakin et 
al. 2009). This model predicts the off-set of the ECAP 
zone toward the opposite direction of the initial defla- 
gration ignition. This second scenario could actually be 
regarded as an extreme case of the first scenario in terms 
of the distribution of the initial deflagration bubbles. 

Although these two scenarios are qualitatively differ- 
ent to one another for the origin of the ECAP zone, the 
resulting configuration (e.g., the off-center ECAP zone) 
may be similar to what we derived in this paper. In 
the second scenario, however, the detonation should pro- 
duce a large amount of 56 Ni in order to produce the 
ECAP zone more massive than ~ 0.1M@ (Meakin et 
al. 2009), at least their model sequence presented to 
date (e.g., - 1.1M Q of 56 Ni in their 2D models). This 
scenario therefore may not be favored as the origin of 
the ECAP zone for the particular case of SN 2003hv, 
which is relatively faint as a normal SN la (Leloudas et 
al. 2009). This second scenario is still a possibility to ex- 
plain the wavelength shift in the nebular emission lines 
seen in bright SNe la, which also show the signature of 
the asymmetry in the ECAP zone (§5 & §6.1). Although 
we expect that the Gravitationally Confined Detonation 
(GCD) model results in bright SNe la if the detonation is 
to produce the ECAP zone, the detailed relation between 
the brightness and the mass of the ECAP zone should be 
dependent on details of the initial conditions such as the 
ignition density (which depends on the accretion rate of 
a white-dwarf before the explosion; e.g., Nomoto 1982; 
Nomoto et al. 1984). Further study based on the GCD 
scenario is therefore important. 

We suggest that future observations of late-time SN 
la emission can provide important information to under- 
stand the details of the explosion mechanism(s). Here- 
after, we discuss importance of the following (future) ob- 
servations, especially in view of the available theoretical 
ideas as mentioned above. 

• Expanding the sample of the line velocity shift mea- 
surement, especially in the optical wavelengths. 

• Studying the detailed line profiles, especially in the 
NIR. 

• Investigating the mass of the ECAP zone. 

6.3.1. Expanding the sample for the line velocity shift 
measurement 

A scries of hydrodynamic simulations based on the 
first scenario provide interesting predictions (Kasen et 
al. 2009). More asymmetric distribution of the initial 
deflagration results in a larger amount of 56 Ni produced 
by the detonation, yielding brighter SNe. For a larger 
degree of the off-set in the initial deflagration bubbles, 
we should also expect a larger degree of the off-set in 
the ECAP zone produced by the deflagration (Maeda et 
al. 2009c). We therefore expect that brighter SNe la 
show a larger off-set velocity in the ECAP zone in this 
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scenario. Since the wavelength shift seen in nebular spec- 
tra of individual SNe la is a result from a combination 
of the degree of the asymmetry and the viewing angle, 
this relation should be observationally investigated in a 
statistic way with a large number of SNe la. We suggest 
to investigate the distribution of the line- velocity (Figure 
10) as a function of the luminosity of SNe la (or so-called 
ATO15 in the light curve shape which correlates with the 
luminosity). The scenario predicts that brighter SNe la 
should show a wider distribution in the velocity space, 
while fainter SNe la should show more narrowly peaked 
distribution at zero velocity. 

The GCD model is an extreme end of the model se- 
quence, predicting bright SNe la. The mechanism to 
create the ECAP zone is different from the first scenario. 
Thus, it may well give the distribution of the velocity 
shift as a function of the luminosity different from the 
first scenario. 

These issues cannot be addressed with the current sam- 
ple of SNe la for which nebular spectra are available - 
the current data are consistent with the idea that SNe 
la have the generic off-set as derived for SN 2003hv (Fig. 
10), but it does not reject a possibility that the off-set 
velocity is dependent on the luminosity because of the 
small sample. Increasing the number of the sample will 
tell us (1) whether brighter SNe la have larger off-set in 
the ECAP zone (a test for the first explosion scenario), 
and (2) whether the distribution of off-set of the ECAP 
zone is explained by a single scenario or hybrid scenarios, 
especially in bright SNe la (a test for the GCD model). 

6.3.2. Detailed Line Profiles especially in the NIR 

In this paper, we have focused on the line wave- 
length shift, but have not tried to fit the details of the 
line profiles. Hoflich et al. (2004) and Motohara et 
al. (2006) suggested to use the detailed line profile of 
[Fe II] 1.644 fj,m to probe the distribution of electron 
capture products (i.e., the ECAP-zone in our model). 
Hoflich et al. (2004) suggested that the apparently flat- 
topped profile of this line in SN 2003du is a signature 
of a "hole" in the region emitting the [Fe II]. They at- 
tributed this hole to the existence of an inner 56 Ni-free, 
electron capture products-rich region, since such a region 
lacks the heating source when the ejecta become trans- 
parent to 7-rays. The same suggestion was also made for 
SN 2003hv based on better data (Motohara et al. 2006). 
There are also other possibilities for the flat-topped pro- 
file, e.g., the infrared catastrophe (Leloudas et al. 2009). 

Understanding the origin of the details of the line pro- 
files requires further studies (e.g., Hoflich et al. 2006). 
For the "hole" interpretation of the flat-topped profiles, 
it should be examined whether such a configuration can 
result from hydrodynamic explosion models. 

An intensive model survey based on the first scenario 
has been presented by Kasen et al. (2009) in 2D. Ac- 
cording to their model sequence, we would expect that 
more asymmetric distribution of the initial deflagration 
bubbles could result in a larger amount of 56 Ni produced 
by the detonation near the zero- velocity. As such, we 
expect that bright SNe la tend to show the peaked line 
profiles. Faint SNe la may tend to show the flat-topped 
profile (sec below) . The most solid case showing the flat- 
top profile in the NIR [Fe II] so far is SN 2003hv. This 
is a relatively faint SN la, thus being consistent with 



this theoretical expectation. Also, the bright SN 1991T 
showed the peaked profile (Bowers et al. 1997) and it is 
also consistent with the expectation. 

Apart from the above prediction, an issue remains 
whether the first scenario can produce the flat-topped 
profiles. The current state-of-art 3D deflagration model 
according to this scenario indicates some amount of un- 
burned C+O is mixed down to the center (e.g., R6pkc 
et al. 2007a). These materials may be later detonated 
to produce 5 %i (Ropke & Niemeyer 2007b; Kasen et al. 
2009). Thus, it is expected that the ECAP zone and 
the HD zone are mixed, rather than form the distinctly 
layered structure. An interesting question is whether fu- 
ture simulations can lead to the " hole" configuration with 
some initial conditions for the thermonuclear bubbles, 
and whether different initial conditions can result in a 
different degree of the mixing (to explain the diversity in 
the detailed profile). 

In the second, GCD scenario, the ECAP zone would 
not suffer from the mixing process since it is created by 
the detonation. It is expected to be always a " hole" , but 
the size of the zone may well depend on the strength of 
the initial deflagration. For the weaker deflagration, the 
density of the white-dwarf at the detonation is higher, 
and thus the size and mass of the ECAP zone are larger 
(Meakin et al. 2009). This raises an interesting possibil- 
ity that the variation in the size of the ECAP zone may 
explain the variation in the detailed line profile; brighter 
SNe la are expected to more likely produce a flat-top 
profile. This tendency is different from the expectation 
from the first scenario. 

In our kinematic model, we have the "hole" in the 
form of the ECAP-zone, which is assumed to be macro- 
scopically separated from the outer 56 Ni-rich HD/LD- 
zones. In this configuration, our model makes the follow- 
ing predictions; (1) Emission lines formed in the ECAP- 
zone have narrow widths (~ 3, 000 km s _1 ), and have a 
peaked profile (with respect to the "off-set" wavelength). 
This is exemplified by [Ni II] A7378, [Ni II] 6.634 /mi, 
[Ni III] 7.350 /im. (2) Emission lines originating from 
the HD-zone have flat-topped profiles. Typical examples 
are [Fe II] A7151, [Fe II] 1.257 /j,m, [Fe II] 1.644 /j,m, and 
[Co III] 11.88 /im. (3) Emission lines dominated by the 
LD-zone have broad, peaked profiles, since the size of the 
"hole" is much smaller than the emitting radius. Typi- 
cal examples are [Fe III] 2.218 /im, and [Fe III] 2.904 /im. 
The [Fe III] blend at 4700A also has the same tendency, 
although the blended nature of this feature makes the 
situation complicated. Some lines show a combination 
of these characteristics; for example, [Fe II] 17.93 /im 
is predicted to show a peaked, relatively broad profile 
with a wavelength shift depending on the viewing angle 
(§5.4). Thus, looking at various emission lines, we will 
be able to obtain different, independent information on 
the explosion physics. 

The best lines for investigating the details of the dis- 
tribution of 56 Ni and the hole (i.e., the ECAP zone) are 
[Fe II] 1.257 /im and [Fe II] 1.644 /im, thanks to their 
low excitation temperature and the relatively isolated 
nature of these lines. Adding to this, [Ni II] A7378 and 
some Mid-IR lines (Table 4) can give the direct view 
of the ECAP zone distribution. Thus, expanding the 
sample of late-time NIR spectra is highly encouraged, 
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preferentially with the optical and (if possible) mid-IR 
data taken almost simultaneously. This will tell us (1) 
whether brighter SNe la have a larger amount of the low- 
velocity 56 Ni (a test for the first explosion scenario), (2) 
whether the fainter SNe tend to show the flat-topped pro- 
file (a constraint on the distribution of the deflagration 
bubbles), and (3) whether the brightest SNe la show the 
hole of the 56 Ni distribution (a test for the GCD model 
as the origin of the brightest SNe la). 

6.3.3. Mass of Neutron-rich Fe-peak Elements 

In §6.2, we have mentioned an importance of multi- 
epoch spectroscopy to clarify details of thermal processes 
within SN la nebulae. This is also important to clarify 
the explosion physics. 

Different scenario should result in the different mass 
of the ECAP zone. The first scenario predicts a smaller 
mass of the ECAP zone for brighter SNe la (since the 
weaker detonation is followed by the stronger detonation; 
Kasen et al. 2009). The trend should be opposite in the 
GCD model sequence (since the neutron-rich materials 
and 56 Ni are both created by the detonation). Deriving 
the mass of neutron-rich materials (i.e., ECAP zone) as a 
function of the SN la peak luminosity will therefore give 
us a strong test for these scenarios. No significant varia- 
tion of the mass of neutron-rich materials as a function of 
the luminosity (or Amis) has been found so far (Mazzali 
et al. 2007). Understanding the thermal properties will 
provide the better accuracy in the mass estimate so that 
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